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Salt  transport  in  the  near-surface  layer  in  the  monsoon-influenced 
Indian  Ocean  using  HYCOM 
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[i]  Salt  advection  in  the  top  5  m  layer  of  the  monsoon- 
influenccd  Indian  Ocean  has  been  studied  with  high  resolu¬ 
tion  HYbrid  Coordinate  Ocean  Model  (HYCOM)  simulated 
sea  surface  salinity  (SSS)  during  2003-06.  The  effect  of 
monsoon  reversals  and  roles  of  Ekman  drift  and  geostrophic 
currents  are  emphasized.  The  Ekman  drift  component  is 
found  to  be  more  important  than  the  geostrophic  component 
in  exchange  of  salt  in  the  basins;  an  indication  of  the  influence 
of  winds  and  the  associated  surface  currents.  The  strongest 
differences  in  the  salt  advection  terms  occurred  in  the  Bay 
of  Bengal  (BoB),  an  indication  of  the  significant  role  of 
advection  for  the  distribution  of  freshwater  and  the  Ekman 
processes.  The  Inter-Tropical  Convergence  Zone  (ITCZ)  is 
seen  to  influence  advection  patterns  via  Ekman  drift  along 
the  equatorial  region.  This  results  in  the  creation  of  tw  o  salt 
advection  bands  whose  positions  vary  seasonally.  This  study 
will  help  to  interpret  data  from  the  recently  launched  Soil 
Moisture  Ocean  Salinity  (SMOS)  and  future  Aquarius  salin¬ 
ity  missions  for  the  Indian  Ocean.  Citation:  Nyadjro,  E.  S., 
B.  Subrahmanyam,  V.  S.  N.  iMurty,  and  J.  F.  Shnver  (2010),  Salt 
transport  in  the  near-surface  layer  in  the  monsoon-influenced 
Indian  Ocean  using  HYCOM,  Geophys.  Res.  Lett.,  37,  LI 5603, 
doi:  1 0.1 029/20 10GL044 127. 

1.  Introduction 

[2]  The  Indian  Ocean  shows  significant  variability  in 
the  distribution  of  sea  surface  salinity  (SSS),  especially  in 
the  two  northern  basins,  the  Bay  of  Bengal  (BoB)  and  the 
Arabian  Sea  (AS)  [Murty  et  al .,  2004;  Joseph  and  Freeland, 
2005].  The  seasonal  reversal  of  monsoon  winds  and  the 
differences  in  precipitation  (P),  evaporation  (E)  and  river 
runoff  (R)  makes  the  BoB  a  net  precipitation  basin  and  the 
AS  a  net  evaporative  basin  [Sengupta  et  al .,  2006].  Also 
the  northern  Indian  Ocean  shows  larger  SSS  variability  than 
the  southern  Indian  Ocean  [Rao  and  Sivaknmar ,  2003].  In 
order  to  balance  salt  and  mass,  there  is  the  need  for  exchange 
between  the  basins.  The  dynamics  of  ocean  salt  distribution 
have  been  studied  in  some  works.  Delcroix  and  Henin 
[1991]  analyzed  and  discussed  equations  for  SSS  changes 
in  the  tropical  Pacific  Ocean.  They  concluded  that  the 
meridional  Ekman  salt  transport  played  a  more  significant 
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role  in  salinity  distribution  than  zonal  salt  transport.  Addi¬ 
tionally  it  was  important  for  the  displacement  of  low  saline 
water  from  regions  where  E  exceeded  P  such  as  the  ITCZ 
and  the  South  Pacific  Convergence  Zone  (SPCZ). 

[3]  Basin-wide  scale  studies  of  SSS  dynamics  in  the 
Indian  Ocean  have  been  a  challenge  and  limiting  as  a  result 
of  paucity  of  observations  of  SSS.  Thus,  the  mechanisms  for 
the  transport  of  salt  in  the  Indian  Ocean  are  not  fully 
understood.  Global  HYCOM  simulated  SSS  data,  validated 
with  in  situ  observations  by  Heffner  et  al.  [2008],  is  used  in 
order  to  make  a  basin-scale  study  of  SSS  dynamics  in  the 
Indian  Ocean. 

[4]  Delcroix  and  Henin  [1991]  expressed  the  time- 
averaged  salt  conservation  equation  for  the  mixed  layer  as: 


(1) 

where  S  is  salinity,  C/thc  zonal  component  of  velocity,  Fthe 
meridional  component  of  velocity,  h  the  depth  of  the  upper 
mixed  layer  (here  only  the  surface  HYCOM  layer,  which 
corresponds  to  the  depth  range  0-5  m,  is  used),  S2=h  the 
salinity  just  below  the  mixed  layer  depth,  H  the  Heaviside 
step  function,  E  the  evaporation,  and  P  the  precipitation.  In 
the  HYCOM  SSS,  the  E-P  and  river  run  off  components  arc 
included.  Vertical  advection  is  usually  not  considered  as  the 
mixed  layer  salinity  is  assumed  to  be  uniform  as  that  of  SSS. 
For  this  reason,  equation  (1)  could  be  considered  valid  for 
SSS  [Delcroix  and  Henin,  1991].  In  furtherance  of  this,  we 
examine  particularly,  the  first  two  terms  of  equation  (1)  and 
their  effects  on  SSS  transport.  Because  of  the  significant 
effect  of  Ekman  forcing  by  monsoonal  winds  and  its  gen¬ 
eration  of  monsoonal  currents  in  the  Indian  Ocean  [Snankar 
et  al.,  1996;  Murtugudde  et  al .,  2007],  we  split  the  hori¬ 
zontal  advection  into  Ekman  drift  and  geostrophic  current 
components  in  order  to  understand  their  relative  roles, 
magnitudes  and  effects  in  the  transport  of  SSS  in  the  Indian 
Ocean.  We  believe  that  the  basin-wide  details  of  such 
studies  are  presently  inexplicit  in  literature. 

[5]  This  paper  hopes  to  contribute  towards  understanding 
the  dynamics  of  salt  distribution  in  the  Indian  Ocean.  The 
strong  annual  variability  in  winds  and  water  circulation,  and 
its  impact  on  the  distribution  of  SSS  also  makes  the  Indian 
Ocean  a  good  laboratory  for  the  calibration  and  validation 
of  the  anticipated  satellite-derived  SSS  from  the  recently 
launched  European  Space  Agency’s  (ESA)  Soil  Moisture 
and  Ocean  Salinity  (SMOS)  and  NASA  and  Space  Agency  of 
Argentina’s  (Comision  Nacional  dc  Actividades  Espacialcs, 
CONAE)  Aquarius/SAC-D  satellite  missions.  The  1/12° 
high  horizontal  resolution  of  the  HYCOM  model  and  its 
smooth  transition  to  a  terrain-following  coordinate  in  coastal 
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Figure  1.  Four-year  mean  HYCOM  salinity  with  boxes  denoting  the  averaging  sub-regions  for  the  study;  Arabian  Sea 
(AS),  Bay  of  Bengal  (BoB),  Sri  Lanka  (SL),  and  Equatorial  (EQ)  regions. 


regions  will  especially  help  in  the  retrieval,  calibration  and 
validation  of  satellite  data  in  coastal  areas  where  sensor 
L-band  measurements  arc  affected  by  nearby  lands  and 
significant  variability  in  river  plumes.  Further,  the  river  run 
off  from  the  major  and  minor  rivers  of  the  Indian  Ocean, 
particularly  into  the  Bay  of  Bengal,  is  included  in  the 
HYCOM  SSS  simulations. 


2.  Data  and  Methods 

[6]  This  study  uses  the  4  year  period  (2003-2006)  monthly 
SSS  from  the  global  HYbrid  Coordinate  Ocean  Model 
(HYCOM)  with  1/12°  horizontal  resolution,  a  description  of 
which  is  provided  and  validated  with  in  situ  observations  by 
Heffner  et  al.  [2008].  The  height  of  the  sea  surface  above  a 
level  surface  C,  was  obtained  from  TOPEX/Poseidon,  and 
Jason- 1  altimetry.  The  surface  wind  stress  r  =  paQo¥  was 
computed  from  the  Florida  State  University’s  (FSU)  Center 
for  Ocean-Atmosphere  Prediction  Studies  (COAPS)  1°  x  1° 
gridded  pseudo-stress  (y/)  data  with  pa  =  1.2041  kg/m3  and 
Co  =  0.0015  [Bourassa  et  al. ,  2005].  Gcostrophic  compo¬ 
nent  of  the  horizontal  salt  transport  was  estimated  from 

^4*  =  ~U8  an<^  Vgs  =  (2) 


where  S  is  salinity,  Ug  and  Vg  arc  surface  gcostrophic  cur¬ 
rents  computed  respectively  from 


<4 


(3) 


where  g  is  acceleration  due  to  gravity  (9.8  ms  2),  /  is  the 
Coriolis  parameter  (20sin 0). 


[7]  Ekman  wind  drift  component  of  the  horizontal  salt 
transport  was  estimated  from 


OS  OS 

Ues  —  Le~^r-  and  Ves  —  \  e 

Ox  ay 


(4) 


where  Ue  and  Ve  are  the  surface  Ekman  drift  computed 
respectively  from  Ue  =  ^  and 


Ve 


(5) 


whereas  ‘h’  is  5  m  depth,  selected  as  the  depth  of  stability 
maximum  in  the  salinity'  stratified  layer,  p0  is  the  water 
density  (1023  kg  m3).  Finally,  surface  salinity  transport  is 
computed  from 


Fs  =  vHAS  (6) 

Where  Fs  is  the  salt  transport  in  psu  m2s_1,  v  is  the  mag¬ 
nitude  of  HYCOM  current  velocity  in  ms-1,  //  is  the  first 
model  layer  depth  in  the  HYCOM  model  (~5  m)  and  AS  is 
the  salinity  anomaly  in  psu.  AS  is  determined  from  AS  - 
S  -  S,  where  S  is  the  annual  mean  SSS. 


3.  Results  and  Discussion 

[8]  Four-year  mean  HYCOM  SSS  (Figure  1)  shows  four 
main  regions  of  salinity  distribution.  The  highest  and  lowest 
values  are  recorded  in  the  AS  and  BoB  respectively.  The 
region  south  of  5°S  have  values  between  those  of  the  AS 
and  BoB.  The  seasonal  difference  in  the  AS  and  BoB  is  well 
simulated  by  the  HYCOM  and  consistent  with  observational 
data  [Rao  and  Sivakumar ,  2003].  The  differences  between 
the  two  basins  can  be  attributed  to  freshwater  forcing 
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Figure  2.  Four-year  mean  variation  of  horizontal  salt 
adveetion  terms  in  (a)  AS,  (b)  BoB,  (e)  SL,  and  (d)  EQ. 


especially  during  the  SW  monsoon  (June-Scptcmber)  when 
precipitation  exceeds  evaporation  in  the  BoB  and  vice  versa 
in  the  AS.  The  low  salinity  tongue  along  the  10°S  latitude  is 
a  consequence  of  the  Indonesian  Throughflow  (ITF)  and  the 
flow  of  BoB  low  salinity  waters  [Sengiipta  et  al. ,  2006]. 

[9]  The  contributions  from  the  salinity  adveetion  terms 
(Ues,  Ugs,  Ves  and  Vgs)  in  the  different  sub-regions  show 
seasonal  and  spatial  variations  (Figures  2a-2d).  The  stron¬ 
gest  differences  in  the  Ekman  component  arc  seen  in  the 
BoB,  where  the  minimum  occurred  in  May  and  November 
(southward  Ekman  component),  and  the  maximum  occurred 
in  August  (eastward  Ekman  component),  emphasizing  the 
influence  of  freshwater  in  salt  transport  especially  during  the 
months  of  monsoon-transition  and  SW  monsoon. 

[10]  The  Ekman  components  of  salt  transport  are  more 
significant  and  4  times  higher  than  the  geostrophie  com¬ 
ponents  (Figure  3).  Summarily,  the  meridional  salt  advec- 
tion  has  a  greater  influence  on  SSS  transport  than  zonal  salt 
adveetion  (Figure  3).  During  the  NE  monsoon  season 
(November-February),  there  is  not  much  difference  in  the 
four  components  of  the  SSS  adveetion  in  the  four  regions 
(Figures  1  and  2)  except  during  November  -  December 
when  salinity  transport  was  dominated  by  the  southward 
Ekman  component  in  the  BoB  (Figure  2b).  This  transports 
low  saline  waters  from  the  BoB,  leading  to  increased  SSS 
along  western  BoB.  During  this  season,  winds  are  pre¬ 
dominantly  from  the  NE,  and  the  Ekman  drift  is  westward 
and  stronger  in  the  western  BoB  [Snankar  et  at.,  1996],  re¬ 
flecting  in  the  adveetion  patterns  during  this  time.  Although 
the  dry  NE  winds  result  in  higher  evaporation  (another 
reason  for  high  SSS  in  western  BoB),  the  eastern  BoB  water 
is  still  low  in  SSS  during  winter  following  the  influx  of 
summer’s  river  runoff.  In  January  and  February,  salinity 
adveetion  in  the  BoB  slows  down,  with  the  geostrophie 
components  being  more  important,  an  indication  of  weak¬ 
ened  wind  effect.  NE  monsoon  salinity  adveetion  is  lower  in 
the  AS  compared  to  the  BoB,  with  the  significant  adveetions 
being  southward  along  the  Somali  current  region  and 
dominated  by  Ekman  drift.  The  rest  of  the  AS  region  shows 
patchy  geostrophie  contributions,  possibly  an  influence  of 
intermittent  eddies.  These  eddy- influenced  transports  also 


occur  in  the  southeastern  Indian  Ocean  where  the  winter 
adveetion  is  predominately  westward  geostrophie  flow. 
Although  / Ian  and  McCreary >  [2001]  observed  Ekman  drift 
in  this  region  to  increase  at  the  peak  of  the  NE  monsoon, 
local  Ekman  pumping  forces  the  geostrophie  flow.  NE 
monsoon  SSS  adveetion  around  the  Sri  Lankan  region  in¬ 
dicates  the  contributions  and  influence  of  different  mecha¬ 
nisms.  At  the  start  of  the  NE  monsoon  season,  the  Ekman 
component  reduces  in  magnitude  and  reverses  from  east¬ 
ward  to  westward,  transporting  water  into  the  AS.  The 
geostrophie  component  becomes  southward  and  lesser  in 
magnitude  than  the  Ekman  component.  At  the  peak  of  the 
NE  monsoon,  all  the  components  are  about  the  same  mag¬ 
nitude  except  the  Ekman  component  which  becomes  very 
significant  (Figure  2).  Geostrophie  contribution  to  SSS 
adveetion  around  the  equatorial  region  was  low  compared  to 
the  Ekman  component  which  shows  two  distinct  continuous 
bands,  with  the  SW  monsoon  adveetion  being  stronger  than 
the  NE  monsoon  (Figures  3e-3h).  During  winter,  there  is  a 
band  of  positive  salt  adveetion  on  either  side  of  the  equator, 
with  the  upper  band  having  a  northward  flow  whiles  the 
lower  band  has  a  southward  flow.  During  the  summer,  both 
bands  lie  on  the  southern  side  of  the  equator  with  the  upper 
band  flowing  northward  while  the  lower  band  flows 
southward,  a  reverse  of  the  winter’s  pattern.  The  pattern 
could  be  linked  to  the  dynamics  of  the  ITCZ.  During  the  NE 
monsoon,  the  ITCZ  migrates  southward,  and  strengthens  the 
NE  trade  winds  that  enhance  evaporation.  The  westward 
extent  of  these  bands  could  also  be  attributed  to  the  pathway 
and  influence  of  the  ITF. 

[11]  The  seasonal  cycle  of  salinity  transport  estimated  in 
the  upper  5  m  layer  from  equation  (6)  using  the  magnitude 
of  current  velocity  and  the  SSS  anomaly  (monthly  HYCOM 
SSS  -  annual  mean  HYCOM  SSS)  in  each  grid  point  is 
presented  in  Figure  4.  We  present  the  bi-monthly  maps  of 
salinity  transport  from  October  (Figures  4a-4f).  This  SSS 
transport  considers  the  model  SSS  (the  E-P  and  river  run  off 
components  are  included)  and  the  model  current  velocity 
and  hence  represents  the  total  magnitude  of  salinity  trans¬ 
port,  without  bifurcating  into  geostrophie  and  Ekman  com¬ 
ponents  of  SSS  transport. 

[12]  In  October,  a  band  of  fresh  water  adveetion  (negative 
salinity'  transport)  is  noticed  along  the  east  coast  of  India  and 
the  eastern  rim  of  the  Bay  and  northern  Sumatra  coast 
(Figure  4a).  This  low  salinity  water  transport  along  the  east 
coast  of  India  is  due  to  the  southward  flowing  East  India 
Coastal  Current  (EICC)  along  the  coast  soon  after  the 
withdrawal  of  the  southwest  monsoon  winds.  McCreary  et  al. 
[1993]  and  Snankar  et  al  [1996]  reported  that  the  southward 
EICC  is  remotely  driven  by  the  downwelling  Kelvin  Waves. 
Interestingly,  during  October,  the  southeastern  Arabian  Sea 
and  south  of  Sri  Lanka  are  the  regions  of  high  salinity 
(positive  salinity  transport).  The  South  Equatorial  Current 
(SEC)  in  the  southern  tropical  Indian  Ocean  transports 
anomalously  high  salinity  waters  towards  west  along  1 2- 1 6°S. 
The  ITF  originating  from  the  Pacific  Ocean  flows  into  the 
SEC.  Potemra  [2005]  mentions  that  the  ITF  transport  is 
maximum  in  general  during  July-October,  the  period  of  SW 
monsoon  and  is  weak  during  NE  monsoon.  Low  salinity 
waters  of  the  Pacific  origin  are  adveeted  into  the  Indian 
Ocean  through  the  ITF.  This  observed  anomalous  high 
salinity  waters  may  be  due  to  the  weakening  of  the  ITF 
during  the  observational  period  (2003-06)  under  the  influ- 
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Figure  3.  Four-year  mean  (top)  East-West  components  and  (bottom)  North-South  components  of  geostrophic  and  Ekman 
drift  of  salinity  advection  during  the  (a,  b,  e,  and  0  Northeast  and  (c,  d,  g,  and  h)  Southwest  monsoon  seasons;  Geostrophic 
and  Ekman  drifts  are  computed  with  the  negative  sign  as  described  in  equations  (2)  and  (4).  Units  are  psu  ms~l. 
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Figure  4.  Bi-monthly  maps  of  salinity  transport  (m2s  *)  with  overlaid  surface  velocity  vectors  (m.s  *)  during  (a)  October, 
(b)  December,  (e)  February,  (d)  April,  (e)  June  and  (f)  August. 


enee  of  the  Indian  Ocean  Dipole  events  [Vinayachandran 
et  al. ,  2007].  Off  northern  Sumatra,  a  pocket  of  low  salin¬ 
ity  water  transport  is  notieed  while  along  the  southern  Sumatra 
a  poeket  of  high  salinity  transport  is  seen.  By  December,  the 
eastern  Bay  transports  low  SSS  waters  towards  Sumatra  eoast 
and  then  westwards  into  the  SEC.  This  observation  is  in 
agreement  with  that  of  Jensen  [2001]  and  Sengupta  et  al. 
[2006]  (Figure  4b). 

[n]  By  February,  the  western  Bay  and  the  east  coast  of 
India  transports  high  salinity  waters  poleward  and  the  North 
Equatorial  Counter  Current  (NECC)  in  the  Southern  Bay 
also  transports  high  SSS  waters  (Figure  4c).  The  South¬ 
eastern  Arabian  Sea  and  Southern  Arabian  Sea  are  now 
replaced  with  the  lower  SSS  and  extends  into  the  south¬ 
western  Arabian  Sea.  A  band  of  low  salinity  water  transport 
is  loeated  at  10°S.  It  is  interesting  to  note  that  both  the 
western  Arabian  Sea  and  western  Bay  of  Bengal  transport 
high  salinity  waters.  However,  the  higher  SSS  in  the  western 
Bay  of  Bengal  is  transported  poleward  and  the  high  SSS 
waters  in  the  western  Arabian  Sea  are  transported  equator- 
ward  by  the  reversing  Somali  eurrent  during  NE  monsoon 
(Figure  4e).  The  open  oeean  regions  (southeastern  Arabian 
Sea  and  south  of  equator)  are  characterized  with  the  low 
SSS  transport.  Part  of  the  low  SSS  transport  in  the  south¬ 
eastern  Arabian  Sea  is  due  to  the  transport  of  the  southward 
EICC  in  December  and  its  further  extension  into  the 
Southeastern  Arabian  Sea  by  February. 

[14]  By  April,  the  low  SSS  ITF  is  seen  as  a  band  along 
10°S  and  its  extension  into  the  Somali  region  along  the  East 
Afriean  Coast  is  evident,  thus  giving  the  evidence  to  the 
contribution  of  ITF  into  the  salt  budget  of  the  Arabian  Sea 


(Figure  4d).  While  the  high  SSS  transport  extends  into  the 
eentral  Arabian  Sea  from  the  Oman  eoast,  the  low  SSS  waters 
of  NE  monsoon  in  the  Southeastern  Arabian  Sea  w  eakened. 
High  SSS  transport  in  the  western  Bay  is  strengthened  by  the 
strength  of  the  poleward  EICC  (Figure  4d).  By  May,  a  band 
of  high  SSS  transport  associated  with  the  eastward  flowing 
Spring  Wyrtki  Jet  is  noticed  in  the  eastern  equatorial  region 
(figure  not  shown).  The  Ekman  transport  of  the  spring 
intermonsoon  winds  brings  in  back  the  low  SSS  waters  from 
the  Southeastern  Arabian  Sea,  into  the  Sri  Lanka  region. 
Compared  to  the  SSS  transport  in  Oetober  (Figure  4a),  the 
SSS  transport  off  Sumatra  is  opposite  in  May  (figure  not 
shown).  With  the  advancement  of  the  southwest  monsoon, 
the  SSS  transport  showed  large  variation  (Figures  4e  and  4f). 

[is]  By  July/August,  the  Southeastern  Arabian  Sea  oc¬ 
cupies  with  high  SSS  waters,  and  Somali  eurrent  transports 
low  SSS  waters  into  the  Arabian  Sea.  These  low  SSS  waters 
in  the  Somali  Currents  are  of  the  ITF  origin  during  April- 
June  (Figures  4d  and  4e).  The  western  Bay  of  Bengal  onee 
again  becomes  the  source  of  low  salinity  waters  from  mon¬ 
soon  precipitation  and  river  runoff.  The  ITF  anomalously 
beeomes  weak  during  July-August  during  the  observation 
period  (2003-06),  as  can  be  seen  from  the  band  of  high  SSS 
transport  along  10°S.  Advection  of  high  SSS  waters  into  the 
Bay  of  Bengal  during  August-September  is  evident  during 
Southwest  monsoon  (Figure  4f). 

4.  Conclusions 

[i6]  The  seasonal  variability  of  SSS  in  the  Indian  Ocean 
as  seen  in  observational  data  is  well  simulated  by  HYCOM; 
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SSS  in  the  Arabian  Sea  is  in  general  higher  than  that  in  the 
Bay  of  Bengal.  Analysis  of  the  contributions  of  Fkman  drift 
and  geostrophic  current  shows  the  former’s  dominance  in 
the  SSS  distribution  in  the  Indian  Ocean  indicating  the 
influence  of  wind  in  the  transport  process.  Although  these 
vary  both  spatially  and  temporally,  on  the  whole  meridional 
transport  dominated  over  zonal  transport.  The  transport  of 
low  saline  waters  into  the  Arabian  Sea  during  the  SW 
monsoon  is  seen  to  originate  from  the  ITF  during  April- 
June.  High  saline  Arabian  Sea  w'ater  is  transported  into  the 
Bay  of  Bengal  along  the  Sri  Lankan  coast  by  strong  geo¬ 
strophic  and  Ekman  contributions.  We  noticed  the  anoma¬ 
lously  weak  ITF  during  the  observation  penod  (2003-06) 
under  the  influence  of  Indian  Ocean  Dipole. 
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